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ABSTRACT 

We  examined  by  simulation  the  use  of  spatial  filtering  for 
enhancing  the  visibility  of  point  targets  and  reducing  the  background 
clutter  in  imagery  obtained  with  a  space-based  staring  mosaic  IR  detec¬ 
tor  array.  An  experimental  digital  image-processing  system  was 
employed  to  interactively  process  Landsat  II  images  containing  simu¬ 
lated  point  targets.  First,  we  obtained  improved  target  detection  and 
clutter  rejection  by  using  filters  that  produced  considerably  more 
high-frequency  enhancement  than  the  conventional  Wiener  filter. 
Second,  convolution  filters  with  large-area  point-spread  functions 
(e.g.,  those  with  11  by  11  elements  or  those  implemented  with  a  2- 
dimensional  FFT  algorithm)  were  more  effective  at  attenuating  the  edges 
of  extended  regions  of  clutter  (e.g.,  those  corresponding  to  clouds, 
bodies  of  water,  etc.)  than  were  those  with  small  areas  (e.g.,  3  by  3 
or  5  by  5  elements). 


RESUME 


Nous  avons  examine  par  simulation  l'aptitude  du  filtrage  spatial 
3  rehausser  la  visibility  des  cibles  ponctuelles  et  A  rSduire 
1'arriSre-plan  dans  des  images  obtenues  par  un  syst&ne  de  surveillance, 
embarqu£  sur  satellites,  employant  une  mosalque  de  d£tecteurs  IR.  Un 
systfeme  experimental  de  traitement  num^rique  d’ images  fut  employe  pour 
traiter  interactivement  des  images  obtenues  par  le  satellite  Landsat 
II  auxquelles  des  cibles  ponctuelles  simul£es  etaient  ajoutees.  Pre- 
miSrement,  les  f litres  qui  accentuent  les  hautes  frequences  sont  prefe¬ 
rables  au  filtre  standard  du  type  Wiener.  Ils  sont  plus  efficaces  pour 
detecter  les  cibles  et  eliminer  l'arridre-plan.  DeuxiSmement ,  les 
filtres  convolutionnels  utllisant  une  fonction  avec  une  grande  surface 
(e.g.  11  x  11  elements  ou  ceux  obtenus  en  utilisant  une  transformation 
de  Fourier  bidimensionnelle)  sont  plus  efficaces  pour  eliminer  les 
contours  des  grands  objets  de  l'arriSre-plan  (e.g.  les  nuages,  les 
lacs,  etc.)  que  les  filtres  avec  une  petite  surface  (e.g.  3  x  3  ou 
5x5  elements). 
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1.0  INTRODUCTION 

When  viewed  from  a  satellite  bearing  a  mosaic  IR  detector  array 
(Refs.  1-5),  small  targets  such  as  aircraft  and  missiles  are  usually 
Imaged  onto  a  single  detector  element.  We  consider  here  the  problem  of 
digitally  enhancing  such  point-target  Images  in  the  presence  of  back¬ 
ground  clutter  (Refs.  6-13),  so  they  can  be  detected  and  tracked  auto¬ 
matically  by  a  space-based  surveillance  system  (Refs.  1-3  and  14).  We 
limit  ourselves  to  the  use  of  spatial  filtering  to  enhance  a  single 
frame  of  imagery.  Temporal  processing  of  a  time  sequence  of  images,  as 
well  as  multispectral  processing,  will  also  be  employed  in  operational 
systems  (Refs.  10,  11  and  13),  but  we  do  not  consider  these  enhancement 
techniques  here. 

The  background  clutter  against  which  point  targets  are  to  be 
enhanced  usually  has  less  high  spatial-frequency  content  than  do  the 
targets.  Therefore,  high-pass  spatial  filtering  can  enhance  the  tar¬ 
gets  relative  to  the  clutter,  and  a  variety  of  such  techniques  have 
been  proposed  (Refs.  10-14).  Some  are  based  on  a  form  of  Wiener  fil¬ 
tering  (Refs.  15-16)  implemented  by  approximating  the  filter  by  a  small 
area  point-spread  function  (PSF)  (usually  one  measuring  3  by  3  ele¬ 
ments)  and  convolving  this  directly  with  the  digital  image. 

First,  whereas  Wiener  filtering  maximizes  the  target-to-clutter 
ratio  in  a  mean-square  error  sense,  it  may  not  be  optimum  for  clutter 
rejection.  Second,  a  filter  with  a  small  point-spread  function  cannot 
efficiently  reject  clutter  that  Is  correlated  over  relatively  large 
areas  of  the  image.  In  particular,  such  a  filter  cannot  effectively 
attenuate  the  edges  of  large  objects  in  a  scene,  such  as  those  corre¬ 
sponding  to  clouds,  lakes  or  rivers.  We  examine  these  two  aspects  of 
point-target  enhancement  and  clutter  rejection  by  digital  computer 
simulation  using  Landsat  satellite  images  containing  artificial  point 
targets. 
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Chapter  2.0  summarizes  the  theory  of  point-target  enhancement  by 
spatial  filtering,  and  compares  direct  spatial-domain  convolution  fil¬ 
tering  with  frequency-domain  filtering  implemented  with  a  two-dimen¬ 
sional  fast  Fourier  transform  (2-D  FFT)  algorithm.  In  Chap.  3.0,  we 
outline  the  digital  image-processing  system  that  we  employed  to  obtain 
the  simulation  results,  and  describe  how  we  generated  the  simulated 
satellite  images.  Simulation  results  comparing  various  forms  of  high- 
pass  spatial  filtering  are  given  in  Chap.  4.0.  In  summary,  the  results 
show  that  more  effective  clutter  rejection  and  point-target  enhancement 
can  be  achieved  by  using  filters  that  (1)  produce  considerably  more 
enhancement  at  high  spatial  frequencies  than  do  Wiener  filters  and  (2) 
have  relatively  large  point-spread  functions  (e.g.,  ones  measuring  11 
by  11  elements  or  2-D  FFT  filters). 

This  work  was  performed  at  DREV  during  the  period  November  1981 
to  July  1982  under  PCN  33A16,  Space-Based  IR  Systems  Study. 

2.0  THEORY 


2.1  Point-Target  Enhancement  by  Spatial  Filtering 

Suppose  we  have  an  Image  that  consists  of  targets  of  interest 
added  to  a  cluttered  background.  If  the  spatial-frequency  contents  of 
the  target  and  background  images  are  different,  we  can  Improve  the 
target-to-clutter  ratio  by  spatial  filtering.  Forms  of  Wiener  fil¬ 
tering  have  been  proposed  as  techniques  for  enhancing  point  targets 
against  cluttered  backgrounds  in  images  obtained  with  a  space-based 
mosaic  IR  detector  array  (Refs.  11-13).  A  Wiener  filter  is  the  linear 
shift-invariant  filter  that  minimizes  the  mean-square  difference 
between  the  intensities  of  the  elements  in  the  restored  image  and  those 
of  the  ideal  signal  of  interest.  In  the  present  situation,  it  will 
optimize  the  amount  of  energy  in  the  image  corresponding  to  the  targets 
relative  to  the  amount  corresponding  to  the  background  clutter. 
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Assume  a  target  image  with  a  power  spectrum  St(u,v)  is  added 
to  a  clutter  image  with  a  power  spectrum  Sc(u,v).  The  result  is 
degraded  by  an  imaging  system  that  introduces  a  blur  characterized  by 
its  Fourier  transform  H(u,v).  In  the  present  application,  H(u,v)  may 
represent  blurring  caused  by  target  or  satellite  motion,  or  by  jitter 
or  other  instabilities  present  in  the  platform  or  the  detector  optics. 
The  variables  u  and  v  represent  the  horizontal  and  vertical  spatial 
frequencies. 

In  the  spatial-frequency  domain,  the  Wiener  filter  (Refs.  15-16) 
is  given  by: 


Q(u,v) 


_ H*(u,v) _ 

jH(u,v)j2  +  Sc(u,v)/St(u,v) 


where  the  asterisk  denotes  the  complex  conjugate.  Assuae  we  have  no 
blurring,  i.e.,  H(u,v)  -  H*(u,v)  «  1  and  that  the  targets  are  points. 
The  latter  means  that  the  power  spectrum  of  the  target  image  is  uni¬ 
form,  i.e.,  St(u,v)  ■  k  where  k  is  a  constant.  The  Wiener  filter 
simplifies  to: 


Q(U’V)  *  k~s"('u,v) 

If  we  further  assume  that  the  point  targets  are  low  enough  in  intensity 
and  few  enough  in  number  that  they  make  a  very  small  contribution  to 
the  total  power  in  the  image  at  all  spatial  frequencies,  i.e., 
Sc(u,v)  »  k,  then 


Q(u,v) 


k 

Sc(u,v) 


The  filter  that  optimizes  the  target-to-clutter  ratio  of  low-intensity 
point  targets  has  a  frequency  response  proportional  to  the  reciprocal 
of  the  power  spectrum  of  the  image. 


UNCLASSIFIED 

4 


The  Wiener  filter  is  an  adaptive  one  that  is  determined  by  the 
power  spectrum  of  the  input  image.  After  its  application,  the  power 
spectrum  of  the  image  will  be  uniform,  i.e.,  it  will  have  been  whit¬ 
ened.  If  a  filter  whitens  the  power  spectrum  of  an  image,  it  can  be 
considered  equivalent  to  the  Wiener  filter  required  for  point-target 
enhancement. 

Wiener  filtering  optimizes  the  targe t-to-clutter  ratio,  but  this 
is  not  necessarily  the  best  criterion  to  use  in  the  present  applica¬ 
tion.  In  particular,  it  may  be  preferable  to  obtain  more  rejection  of 
low-frequency  clutter  at  the  expense  of  an  overall  decrease  in  the 
targe t-to-clutter  ratio.  The  low-frequency  clutter  may  be  difficult  to 
remove  by  other  forms  of  processing,  and  it  could  significantly  hinder 
automatic  detection  and  tracking  of  the  targets.  For  example,  inaccu¬ 
racies  in  registering  successive  image  frames  obtained  at  different 
times  may  prevent  frame  subtraction  (Refs.  17-18)  from  cancelling  the 
edges  of  extended  areas  of  clutter.  Retention  of  such  edges  can  pres¬ 
ent  particular  difficulties  for  certain  types  of  tracking  algorithms 
(Ref 8.  19-20),  and  as  a  result  the  edges  may  be  classified  as  false 

target  tracks. 

In  selecting  filters  to  enhance  high  frequencies  or  attenuate 
low  ones,  it  is  important  to  ensure  that  the  frequency  response  of  the 
filter  is  as  smooth  and  gradual  as  possible.  High-pass  filters  with 
sharp  low-frequency  cutoffs  can  introduce  artifacts,  such  as  ringing  at 
the  edges  of  objects  or  regions.  One  type  of  high-pass  filter  that 
satisfies  this  smooth-cutoff  criterion  and  which  permits  an  easy 
adjustment  of  the  amount  of  high-frequency  enhancement  is  the  circu¬ 
larly  symmetric  class  of  filters  defined  by: 


F(ui)  -  un 
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where  io2  -  u2  +  v2  and  n  is  a  real  number.  For  certain  values  of  n 
(typically  between  2  and  3)  these  filters  approximately  whiten  the 
power  spectra  of  the  landsat  II  Images  described  in  Chap.  3.0  and, 
therefore,  represent  approximations  to  the  point-target  Wiener  filter. 
We  examine  experimentally  this  class  of  filters  in  Chap.  4.0. 

A  variety  of  processing  stages  will  be  employed  in  a  space-based 
IR  surveillance  system.  In  addition  to  the  spatial  filtering  consid¬ 
ered  here,  they  include  temporal  processing,  multispectral  processing, 
and  target  detection  and  tracking.  The  optimum  spatial  filtering  can¬ 
not  be  determined  without  also  including  the  other  forms  of  processing. 
However,  we  will  show  in  Chap.  4.0  that  in  optimizing  the  overall  sys¬ 
tem  performance,  one  should  consider  spatial  filters  which  are  quite 
different  from  the  Wiener  filter  that  maximizes  the  target-to-clutter 
ratio. 

2.2  Frequency-Domain  and  Spatial-Domain  Filtering 

The  spatial  filters  described  in  Sect.  2.1  were  defined  in  the 
spatial-frequency  domain.  They  can  be  implemented  exactly  by  calculat¬ 
ing  the  2-D  FFT  of  an  image,  multiplying  by  the  specified  filter  fre¬ 
quency  response,  and  then  calculating  the  inverse  2-D  FFT. 

Such  filters  can  also  be  implemented  by  direct  spatial-domain 
convolution.  In  this  case,  the  filter  is  defined  by  its  PSF,  which  is 
the  result  obtained  when  the  frequency-domain  filter  is  applied  to  an 
image  containing  a  single  point  against  a  zero  background.  The 
resulting  PSF  of  the  frequency-domain  filter  is  then  truncated  to  a 
convenient  size  and  convolved  with  the  image.  However,  as  we  examine 
experimentally  in  Sect.  4.3,  the  truncation  of  the  filter  PSF  results 
in  a  frequency  response  that  only  approximates  the  desired  one,  and 
this  can  introduce  undesired  artifacts  or  deteriorate  the  filter  per¬ 
formance. 


I 


UNCLASSIFIED 

6 


3.0  DIGITAL  IMAGE-PROCESSING  SIMULATION 

/ 

3.1  Experimental  System  and  Process!*.?  Techniques 

The  results  described  in  Chap.  4.0  were  obtained  by  simulation 
using  the  experimental  digital  image-processing  system  at  DREV 
(Ref.  21).  A  variety  of  general-purpose  image-processing  algorithms 
are  available,  and  these  are  controlled  interactively  by  an  operator 
who  gives  commands  on  keyboards  or  with  a  track  ball.  The  main  proc¬ 
essing  operations  used  to  obtain  the  present  results  Include  those  for 
performing  2-D  FFT  filtering,  2-D  spatial-domain  convolution,  intensity 
transformation,  power-spectrum  analysis,  display  enhancement  and 
graphics  generation. 

Figure  1  shows  a  block  diagram  of  the  system.  It  is  based  on  a 
multi-user  PDP  11/70  computer,  with  standard  peripherals,  which  con¬ 
trols  two  high-speed  processing  elements.  One  of  the  processing  ele¬ 
ments  is  a  Floating  Point  Systems  Inc.  AP120B  programmable  array  trans¬ 
form  processor,  and  the  other  one  is  a  Comtal  Vision  0ne/20  interactive 
image-processing  and  display  system. 

The  array  processor  operates  at  12  million  operations  per  second 
(38-bit  words)  and  permits  a  variety  of  processing  algorithms,  such  as 
2-D  FFT  filtering,  to  be  readily  programmed  in  FORTRAN  language  and 
applied  in  several  seconds  to  256x256-eleraent  images. 

The  Comtal  system  operates  at  141  million  operations  per  second 
(8-bit  words)  and  permits  images  to  be  digitized  and  processed  at  real¬ 
time  TV  rates.  Live  TV  video  is  digitized  with  a  spatial  resolution  of 
512  by  512  elements,  and  an  intensity  resolution  of  256  levels,  at  a 
rate  of  30  frames  per  second.  Internal  pipeline  processors  permit 
convolution  (i.e.,  spatial  filtering  with  a  3x3-element  PSF),  image 


UNCLASSIFIED 

7 


DR EU  TRACE-PROCESSING  SVSTEM 


12  DISKS! 

160  MB 


2  MflG 
TAPES 


DISK 
1.2  MB! 


LINE 

PRINTR 


P 

D 

P 

1 

i 

/ 

7 

0 

C 

P 

U 


USER 

USER 

TTYS 

\  ? 

TTY 

ARRAY 

PROC. 

u 

GPIOP 

H 

DISK 
300  MB 


TU 

CAMERA 

COLOR 

W 

[ 

IMAGE 

FILM 

PROC. 

WRITER 

.  1  . 

kVBRD, 

VTR 

TBALL 

EQPMT. 

STANDARD 

PERIPHERALS 


I MAGE-PROC ESS  I NG  DEU I CES 


FIGURE  1  -  Block  diagram  of  the  experimental  digital  image-processing 
system  used  to  investigate  the  spatial  filters 


combination  (i.e.,  addition,  subtraction,  multiplication  or  division 
of  two  stored  or  live  TV  images),  intensity  transformation,  and  pseudo¬ 
coloring  to  be  implemented  at  the  real-time  TV  rate. 

Images  can  be  rapidly  transferred  between  the  array  processor, 
Comtal  system,  disk  memory  and  magnetic  tape.  They  can  also  be 
recorded  on  35-mm,  4-  by  5-in  Polaroid  and  8-  by  10-In  Polaroid  films 
with  a  Matrix  Instruments  model  4007  color  graphic  camera  interfaced  to 
the  Comtal  system. 
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The  six  original  Landsat  satellite  Images  employed  In  Chap.  4.0 
were  recorded  on  70-mm  positive  transparency  film.  We  placed  each  film 
Image  on  a  light  table  and  viewed  It  with  a  hlgh-resolutlon  vidicon 
camera  made  by  Spatial  Data  Inc.  Four  TV  frames  were  digitized,  aver¬ 
aged  and  stored  In  an  Image  memory.  The  camera  lens  was  covered  and  a 
16-frame  average  was  similarly  acquired  and  stored.  We  then  subtracted 
the  latter  'black  level'  digital  Image  from  the  digitized  film  Image, 
and  minified  this  by  a  factor  of  two  by  replacing  2-  by  2-element 
squares  by  their  averages  to  produce  a  256-  by  256-element  Image. 
Point  targets  were  added  to  this  Image  as  described  In  Sect.  3.2,  and 
the  result  was  processed  with  38-blt  accuracy  In  the  array  processor  by 
using  the  filtering  techniques  described  In  Chap.  4.0.  The  processed 
Images  were  normalized,  displayed  on  the  Comtal  system,  and  recorded  on 
Polaroid  film  with  the  Matrix  camera.  Most  of  the  figures  Included  In 
Chap.  4.0  were  obtained  by  using  a  multiple-image  format  In  which  six 
Images  were  automatically  recorded  on  a  single  sheet  of  8-  by  10-in 
Polaroid  film. 

3.2  Simulation  Approximations 

Imagery  obtained  In  the  far  IR  (8-13  urn)  with  a  satellite-based 
mosaic  detector  which  shows  point  targets  against  a  cluttered  back¬ 
ground  Is  not  currently  available  at  DREV.  We  generated  simulated 
Images  to  permit  us  to  Investigate  the  processing  techniques  that  may 
be  employed  by  a  space-based  IR  surveillance  system.  The  background 
clutter  Is  represented  in  the  simulation  by  Landsat  II  satellite 
imagery  obtained  In  the  near  IR  (0. 8-1.1  pm),  and  we  added  point 
targets  to  selected  Image  elements. 

Objects  or  regions  In  a  scene  will  not  necessarily  have  the  same 
relative  Intensities  when  viewed  with  near-  and  far-IR  radiation.  This 
is  evident  In  Figs.  2(A)  and  2(B)  which  give  images  of  the  city  of 
Guelph,  Ontario,  obtained  simultaneously  In  visible  light  and  In  the 
far  IR  with  a  multispectral  scanner  (MSS)  located  on  board  an  air- 
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FIGURE  2  -  Images  from  an  airborne  nultispectral  scanner  obtained 
simultaneously  with  visible-light  (A)  and  IR  radiation  (B). 
The  power  spectra  (C)  and  (D)  of  the  two  images  are  simi¬ 
lar.  The  inserts  in  the  power  spectra  show  graphs  of  their 
radial  power  distributions. 
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craft.  However,  the  visible  aod  far-IR  Images  have  approximately  the 
same  spatial-frequency  content  as  their  power  spectra  (Figs.  2(C)  and 
2(D))  show.  Since  the  present  processing  Is  based  on  spatial-frequency 
content,  not  directly  on  the  relative  contrast  of  the  different 
regions,  the  near-IR  Images  that  we  use  to  represent  the  far-IR  back¬ 
ground  clutter  should  be  a  reasonable  approximation  for  the  present 
purposes. 

We  obtained  the  power  spectra  shown  In  Figs.  2(C)  and  2(D)  by 
calculating  the  modulus  squared  of  the  2-D  FFT  of  each  image.  In  each 
power-spectrum  Image,  the  zero  spatial  frequency  Is  at  the  center  and 
the  Nyqulst  frequencies  are  at  the  borders.  The  power  spectra  were 
displayed  with  a  4-decade  logarithmic  intensity  scale,  therefore  full 
black  corresponds  to  a  level  80  dB  below  the  maximum.  Since  the  origi¬ 
nal  images  are  real,  the  lower  and  upper  halves  of  the  power  spectra 
are  mirror  Images.  In  the  lower  right  quadrant  of  each  power  spectrum 
we  have  graphed  Its  radial  power  distribution.  This  was  obtained  by 
Integrating  the  power  within  flxed-thlckness  rings  centered  about  zero 
spatial  frequency,  and  plotting  this  as  a  linear  function  of  the  radius 
of  the  center  of  the  ring. 

We  simulated  the  point  targets  by  adding  constants  to  a  number 
of  nonadjacent  Image  elements.  This  represents  targets  that  have  a 
higher  apparent  temperature  than  the  background  and  which  are  much 
smaller  than  the  field  of  view  of  an  Individual  detector  element.  The 
target  images  are  assumed  to  be  small  enough  that  they  always  appear  on 
a  single  detector  element,  and  never  overlap  two  neighboring  ones.  In 
practice,  the  point  targets  will  be  broadened  by  the  finite  resolution 
of  the  imaging  system.  As  will  be  described  In  Sect.  3.3,  we  use 
Images  that  have  been  digitized  with  a  relatively  low  spatial  resolu¬ 
tion  (200  m  on  the  ground).  In  practice,  It  should  be  possible  to 
employ  optics  with  a  much  higher  resolution  than  this,  so  there  would 
be  no  significant  broadening  of  the  point  targets.  However,  depending 
on  the  size  of  the  gaps  (i.e.,  the  Inactive  areas)  between  the  elements 
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of  the  mosaic  IR  detector  array  at  this  spatial  resolution,  large  tar¬ 
gets  (such  as  aircraft)  may  occupy  more  than  a  single  element  In  the 
digital  Image. 

Mosaic  IR  detector  arrays  exhibit  nonuniformities  In  overall 
sensitivity  to  radiation  and  In  dark  current  from  one  element  to 
another.  Processing  techniques  (Refs.  22-23)  can  be  used  to  reduce 
these  effects.  Such  point-to-point  nonuniformities  do  not  occur  In 
Landsat  II  MSS  Images.  It  is  not  known  how  such  nonuniformities  will 
affect  the  ability  of  high-pass  spatial  filtering  to  reduce  clutter  and 
enhance  the  point  targets.  Approximation  of  a  mosaic  detector  array  by 
the  Landsat  MSS  Is  probably  the  most  questionable  assixnption  involved 
in  generating  the  present  simulated  Images. 

3.3  Production  of  the  Simulated  Images 

We  used  images  of  six  different  scenes  acquired  In  Band  7  of  the 
MSS  located  on  board  the  Landsat  II  satellite.  The  images  were  formed 
with  near-IR  radiation  In  the  range  0. 8-1.1  urn,  and  were  recorded  on 
70-mm  black/whlte  positive  transparency  film.  We  selected  1.5-  by  1.5- 
cm  regions  of  each  film  image  and  digitized  these  with  a  spatial  reso¬ 
lution  of  256  by  256  elements  using  the  procedure  described  in 
Sect.  3.1.  Each  element  in  the  resulting  digital  images  corresponds  to 
a  square  on  the  ground  measuring  200  m  on  a  side. 

Figure  3  shows  the  six  digital  images  after  the  point  targets 
have  been  added  as  described  below.  We  used  homomorphic  filtering 
(Ref.  24)  to  improve  the  visibility  of  low-contrast  details  in  the 
display  and  to  compress  the  dynamic  range.  The  resulting  enhanced 
images  are  given  in  Fig.  4.  The  six  images  represent  different  types 
of  high-clutter  backgrounds.  Images  A  and  B  show  an  area  in  the  vicin¬ 
ity  of  Quebec  City  in  summer,  with  and  without  cloud  cover,  wher*-  is 
Images  C  and  D  show  the  Montreal  area  in  winter  and  in  simmer.  A 


FIGURE  3  -  Six  near-IR  images  obtained  with  the  Landsat  II  satellite. 

Each  represents  a  different  type  of  high-clutter  back¬ 
ground,  and  contains  64  simulated  low-intensity  point  tar¬ 
gets.  The  images  are  composed  of  arrays  of  256  by  256 
elements  and  have  a  ground  resolution  of  200  m/element. 


UNCLASSIFIED 

13 


FIGURE  4  -  The  Landsat  II  Images  of  Fig.  3  after  homomorphic  filtering 
to  enhance  low-contrast  information  and  compress  the 
dynamic  range 
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TABLE  1 


CHhPhi  TEFIST  ICS  Of  THE  SIX  LAMOSAT  II  IMAGES 


IMAGE 

LAND SAT  MO 

0P61T 

FRAME 

[ATE 

SUBJECT 

A 

2237814495 

14 

27 

27-07-81 

QUEBEC  Cl  TV  -  CLOUDY 

B 

2239614494 

14 

27 

14-08-81 

QUEBEC  CITY  -  CLEAR 

C 

2212815045 

16 

28 

19-11-80 

MUM T PE PL  -  UINTEP 

D 

2236215021 

16 

28 

11-07-81 

MOHTPEAL  -  SUMMEP 

E 

2236715290 

21 

24 

16-07-81 

EASTMAN.  JAMES  B 0V 

F 

2234915320 

21 

30 

28-0*.-8l 

E'ETPOIT  CITY 

FIGURE  5  -  The  simulated  point  targets  that  were  added  to  each  of  the 
Landsat  II  images  of  Fig.  3.  The  point  targets  have  been 
enlarged  to  an  area  of  3  by  3  elements  to  improve  their 
visibility. 
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region  near  James  Bay  viewed  In  the  summer  is  given  in  Image  E,  and  the 
Detroit  City  area  is  shown  in  Image  F.  Table  I  gives  more  detailed 
information  on  the  six  Landsat  II  images. 

Each  of  the  six  images  shown  in  Figs.  3  and  4  contains  64  point 
targets.  As  illustrated  in  Fig.  5,  the  point  targets  are  arranged  in 
eight  rows  of  eight  targets  each.  The  intensities  of  the  targets 
Increase  from  left  to  right  in  rows  1,  3,  5,  etc.  (counting  the  rows 
from  the  bottom),  whereas  in  rows  2,  4,  6,  etc.  they  increase  in  inten¬ 
sity  from  right  to  left.  In  each  row,  the  target  intensity  varies  in 
seven  linear  steps  from  a  minimum  of  0.024M  to  a  maximum  of  0. 2M,  where 
M  is  the  maximum  intensity  in  the  image.  The  RMS  contrast  of  the  aver¬ 
age  image  is  0.146M,  so  we  have  generated  targets  with  signal-to-clut- 
ter  ratios  between  0.164  and  1.37.  Note  that  in  Figs.  3  and  4  even 
the  most  intense  point  targets  are  almost  invisible  in  the  high-clutter 
areas  of  the  images,  and  are  only  marginally  visible  against  the  bodies 
of  water  that  have  more  uniform  intensity  distributions. 

4.0  ENHANCEMENT  RESULTS 

4.1  Frequency-Domain  Filtering 

We  selected  Landsat  Image  A,  shown  in  Fig.  3(A),  and  applied  to 
it_six  different  high-pass  spatial  filters.  The  filters  were  imple¬ 
mented  in  the  spatial-frequency  domain  with  a  2-D  FFT  algorithm,  and 
had  the  power-law  response  F(w)  ■  u>n  described  in  Sect.  2.1.  Figure 
6(A)  shows  the  original  image  whereas  Figs.  6(B)  to  6(F)  give  the 

results  obtained  by  using  filters  with  frequency  responses  o,2,  u>3 ,  w4 , 
o>6  and  a,3.  Before  displaying  each  processed  image,  we  normalized  its 

contrast  by  setting  the  black  level  of  the  display  to  1^  -  1.5a  and 

its  full-white  level  to  1^  +  6o ,  where  IA  and  a  are  the  mean  and 

the  standard  deviations  of  the  image  intensity  distribution. 
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FIGURE  6  -  Landsat  Image  A  after  five  different  types  of  high-pass 
spatial  filtering  implemented  in  the  spatial-frequency 
domain  with  a  2-D  FFT.  Increasing  the  amount  of  high-fre¬ 
quency  enhancement  improves  the  visibility  of  the  point 
targets  and  gives  better  rejection  of  the  background  clut¬ 
ter. 
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Figure  6  clearly  shews  that  filters  which  produce  a  small  or  a 
moderate  amount  of  high-frequency  enhancement  (e.g.,  id2  oru3)  are  less 
effective  at  reducing  the  background  clutter  and  enhancing  the  point 
targets  than  are  those  that  produce  a  larger  amount  of  it  (e.g. ,  id6  or 
id8).  For  example,  the  u>2  and  <d3  filters  retain  clutter  corresponding 
to  the  edges  of  clouds  and  the  land-water  boundaries.  These  edges  are 
almost  completely  removed  by  the  <d8  and  id8  filters.  Note  the  strong 
enhancement  of  the  point  targets  relative  to  the  remaining  background 
clutter  that  the  latter  two  filters  have  achieved. 

The  improvement  in  target  detectability  is  shown  graphically  in 
Fig.  7.  In  this  figure,  we  plotted  the  number  of  targets  detected  in 
Landsat  Image  A  as  a  function  of  the  frequency  response  of  the  high-pass 
filter  for  four  different  false-alarm  rates.  We  obtained  the  upper 
curve  in  Fig.  7  by  applying  a  threshold  level  to  each  filtered  image, 
adjusting  it  to  obtain  16  false  alarms  in  the  image,  and  then  counting 
the  nuaber  of  point  targets  that  were  correctly  detected.  The  lower 
three  curves  were  similarly  obtained  by  using  detection  thresholds  that 
produced  eight,  four  and  two  false  alarms. 

Increasing  the  amount  of  high-frequency  enhancement  from  id2  to  <d6 
improves  the  target-detection  probability.  Note  that  the  detection 
probability  is  a  nonlinear  function  of  the  ntsnber  of  detected  targets 
since  the  ones  present  in  Landsat  Image  A  have  varying  intensities.  For 
example,  for  the  curve  representing  two  false  alarms,  increasing  the 
number  of  detected  targets  from  35  (for  the  id2  filter)  to  45  (for  the  id6 
filter)  represents  a  substantial  reduction  in  the  minimum  intensity  a 
target  must  have  to  be  detected  with  reasonable  probability. 

Figure  8  gives  the  pewer  spectra  and  the  corresponding  radial 
distribution  functions  of  the  original  and  the  five  high-pass  filtered 
Images  of  Fig.  6.  We  obtained  the  spectra  and  the  distribution  func¬ 
tions  as  described  in  Sect.  3.2.  As  discussed  in  Sect.  2.1,  the  filter 
that  whitens  the  power  spectrum  of  an  image,  i.e.,  that  modifies  it  to 
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FILTER  FREQUENCY  RESPONSE 


FIGURE  7  -  The  target-detection  performances  of  the  five  frequency- 
domain  filters  shown  in  Fig.  6.  We  generated  each  of  the 
four  curves  by  setting  a  single-intensity  threshold  on  the 
filtered  images  so  as  to  obtain  the  specified  number  of 
false  alarms,  and  then  counting  how  many  of  the  64  true 
point  targets  were  detected. 
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FIGURE  8  -  The  power  spectrum  of  the  original  Landsat  Image  A  (A)  is 
shown  in  (B)  to  (F)  after  the  application  of  five  different 
f req uency -domain  high-pass  filters.  The  inserts  in  the  power 
spectra  show  graphs  of  their  radial  power  distributions. 
The  power  spectrum  of  this  image  is  approximately  'whitened' 
by  a  filter  with  a  response  proportional  to  the  square  or 
the  cube  of  the  spatial  frequency. 
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have  equal  power  at  all  spatial  frequencies,  represents  the  Wiener  fil¬ 
ter  required  to  optimize  the  signal- to-noise  ratio  of  low-intensity 
point  targets.  Figure  8  illustrates  that,  except  at  very  low  spatial 
frequencies,  a  filter  with  a  frequency  response  between  <d2  and  <d3 
approximates  the  point-target  Wiener  filter  for  Image  A.  Since  a  much 
improved  target  enhancement  and  clutter  rejection  is  obtained  by  using 
considerably  more  high-frequency  enhancement  (e.g.,  F(<u)  »  <d6),  the 
Wiener  filter  does  not  appear  to  be  optimum  for  the  present  situation. 

Figure  9  shows  the  results  of  applying  the  w6  filter  to  the  six 
Landsat  images  of  Fig.  3.  In  all  cases  the  high-pass  filtering 
strongly  enhanced  the  point  targets  and  suppressed  the  background.  The 
display  normalization  was  the  same  as  that  used  in  Fig.  6. 

In  Fig.  10,  we  attempted  to  extract  the  point  targets  from  the 
background  by  setting  a  single-intensity  threshold  on  the  filtered 
images  of  Hg.  9.  If  an  element  had  an  intensity  belcw  IA  +  6o  it 
was  set  to  black,  and  if  it  equaled  or  exceeded  this  threshold  it  was 
set  to  full  white.  To  improve  the  visibility  of  the  point  targets 
extracted  in  Fig.  10,  we  enlarged  each  target  to  a  3-  by  3-element 
square.  Only  11  false  alarms  are  present  in  the  six  thresholded 
images,  and  266  out  of  the  384  points  were  correctly  detected.  All 
targets  with  an  intensity  of  0.656a  or  higher  were  detected.  Those 
with  an  intensity  of  0.492a  were  detected  with  a  50X  probability,  but 
only  2 X  of  the  0.328a  targets  were  extracted. 

Increasing  the  high-frequency  content  of  the  background  clutter 
made  it  more  difficult  to  detect  the  point  targets.  For  example,  the 
targets  that  were  present  against  the  lower  frequency  backgrounds  shown 
in  Images  B  and  E  were  detected  with  a  higher  probability  than  those 
present  against  the  higher  frequency  backgrounds  of  Images  A  and  F. 
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FIGURE  9  -  The  six  Landsat  images  of  Fig.  3  after  high-pass  filtering 
to  enhance  the  point  targets  and  reduce  the  background  clut¬ 
ter.  The  filter  was  the  same  for  all  cases  and  had  a 
response  proportional  to  the  sixth  power  of  the  spatial 
frequency. 
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FIGURE  10  -  The  six  high-pass  filtered  Landsat  images  of  Fig.  9  after 
setting  a  single-intensity  threshold  at  IA  +  fo  to  detect 
the  point  targets 
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Histograms  showing  the  intensity  distributions  of  the  six  origi¬ 
nal  Landsat  images  of  Fig.  3  are  given  in  Fig.  11.  Most  distributions 
are  bimodal,  with  the  sharp  lew-intensity  peak  corresponding  to  the 
darker  water  areas  and  the  broader  mid-intensity  peak  corresponding  to 
the  land  areas.  As  shown  in  Fig.  12,  after  high-pass  filtering  with 
F(ui)  -  u)6,  all  images  have  an  intensity  distribution  that  is  approxi¬ 
mately  Gaussian. 

This  is  opportune  because  it  permits  us  to  use  normal  statistics 
in  analyzing  and  implementing  the  target-extraction  operation.  In 
particular,  to  detect  the  targets  in  each  of  the  present  high-pass 
filtered  images,  it  was  necessary  to  calculate  only  two  parameters  of 
the  intensity  distribution:  its  mean  and  its  variance.  By  setting  a 
fixed-intensity  threshold  at  IA  +  6o ,  we  obtained  an  acceptable  bal¬ 
ance  between  miss  and  false-alarm  rates  for  all  images,  despite  the 
significant  differences  in  the  characteristics  of  the  background  clut¬ 
ter. 


4.2  Spatial-Domain  Filtering 

As  described  in  Sect.  2.2,  the  implementation  in  the  spatial- 
frequency  domain  of  a  filter  is  equivalent  to  the  convolution  of  its 
PSF  with  the  original  image.  In  this  section  we  calculate  the  PSF  of 
the  <j)6  filter  employed  in  Sect.  4.1,  and  examine  the  effect  of  trun¬ 
cating  it.  In  particular,  we  apply  to  Landsat  Image  A  3x3-,  5x5-, 
7x7-,  9x9-  and  llxll-element  approximations  of  the  <*)6  filter. 

We  generated  the  PSF  of  the  id6  filter  by  calculating  the  2-D  FFT 
of  an  image  containing  a  point  against  a  zero  background,  multiplying 
this  by  the  function  u6  and  performing  the  inverse  2-D  FFT.  The 
resulting  PSF,  enlarged  by  a  factor  of  4,  is  shown  in  Fig.  13.  Only 
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FIGURE  11  -  Histograms  showing  the  distribution  of  intensities  in 
Landsat  Images  A  to  F.  The  sharp  peaks  at  low  intensity 
correspond  to  water  and  the  broader  central  peaks  corre¬ 
spond  to  land. 
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FIGURE  12  -  Histograms  showing  the  distribution  of  intensities  in 
Landsat  Images  A  to  F  after  application  of  a  high-pass 
filter  with  a  response  proportional  to  the  sixth  power  of 
the  spatial  frequency.  The  filtering  has  removed  the 
bimodal  character  evident  in  Fig.  11  and  has  produced 
almost  normal  intensity  distributions. 
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FIGURE  13  -  The  PSF  of  a  frequency-domain  filter  with  a  response  pro¬ 
portional  to  the  sixth  power  of  the  spatial  frequency. 
Each  filter  coefficient  is  displayed  as  a  square  element 
with  an  intensity  proportional  to  its  value.  Zero  appears 
with  a  mid-gray  intensity;  positive  and  negative  coeffi¬ 
cients  have  higher  and  lower  intensities. 


FIGURE  14  -  The  PSF  of  Fig.  13  after  truncation  to  form  five  different 
small-area  approximations 
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the  central  64x64-element  region  of  the  PSF  is  displayed,  and  the  value 
of  each  filter  coefficient  is  represented  by  the  intensity  of  a  4x4- 
eleinent  square.  The  display  has  been  biased  so  that  zero  intensity 
appears  gray  and  positive  and  negative  filter  coefficients  appear, 
respectively,  lighter  and  darker  than  this.  We  used  a  nonlinear  inten¬ 
sity  transformation  to  improve  the  visibility  of  low-intensity  positive 
and  negative  filter  coefficients.  The  filter  is  circularly  symmetric, 
so  it  can  be  represented  by  the  coefficients  within  a  single  quadrant. 
A  listing  of  the  relative  values  of  the  coefficients  in  the  lower  right 
quadrant  of  the  PSF  is  given  in  Table  II.  Note  that  the  positive  and 
negative  filter  coefficients  are  interleaved. 

We  generated  five  filters  with  PSFs  approximating  that  of  the 
one.  These  were  obtained  by  retaining  the  filter  coefficients  within 
3x3-,  5x5-,  7x7-,  9x9-  and  llxl 1-element  regions  of  the  center  of  the 
PSF  of  the  to6  filter,  and  setting  to  zero  all  the  other  coefficients. 
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The  central  coefficient  of  each  PSF  was  then  adjusted  to  obtain  a 
minimum  filter  response  at  zero  spatial  frequency.  Figure  14  shows  the 
five  resulting  truncated  PSFs. 

Figure  15  shows  the  result  of  convolving  the  five  PSFs  shown  In 
Fig.  14  with  Landsat  Image  A,  along  with  that  obtained  by  Implementing 
the  a >6  filter  with  a  2-D  FFT.  We  normalized  the  display  of  each  fil¬ 
tered  image  so  that  full  black  corresponded  to  intensity  I ^  -  2a  and 
full  white  to  intensity  IA  +  6o .  Since  the  total  energy  correspond¬ 
ing  to  the  targets  is  small  compared  to  that  of  the  clutter,  the  nor¬ 
malization  has  the  effect  of  displaying  the  clutter  with  approximately 
the  same  overall  contrast  in  each  image.  The  six  filtered  images  were 
thresholded  at  a  level  of  1^  +  5.5a,  as  described  in  Sect.  4.1,  to 
produce  the  detected  targets  shown  in  Fig.  16. 

The  clutter  remaining  after  the  application  of  a  small-area 
filter  (e.g.,  one  with  3  by  3  or  5  by  5  elements)  Is  more  strongly 
correlated  with  extended  features  and  edges  in  the  original  image  than 
Is  that  remaining  after  the  use  of  a  large-area  one  (e.g.,  the  11x11- 
element  or  the  2-D  FFT  filter).  For  example,  the  edges  of  the  river 
are  clearly  apparent  after  the  application  of  the  3x3-element  filter, 
whereas  these  features  are  barely  visible  after  the  application  of  the 
llxll-element  and  the  2-D  FFT  filters.  As  discussed  In  Sect.  2.1, 
clutter  that  is  correlated  with  the  edges  of  extended  objects  in  a 
scene  may  be  detrimental  to  certain  automatic  target-detection  and 
tracking  algorithms. 

4.3  Comparison  of  the  Spatial  Filters 

We  have  graphed  in  Fig.  17  the  number  of  targets  detected  in 
Landsat  Image  A  as  a  function  of  the  type  of  high-pass  filter  that  was 
used.  Results  are  shown  for  the  five  convolution  filters  considered  in 
Sect.  4.2  as  well  as  for  the  2-D  FFT  filter  employed  in  Sect.  4.1.  The 
graphs  were  obtained  as  described  in  Sect.  4.1  for  fixed  false-alarm 
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FIGURE  15  -  Landsat  Image  A  after  high-pass  filtering  with  five  dif¬ 
ferent  8 pa tia 1-domain  convolution  filters  and  a  frequency- 
domain  one  with  a  response  proportional  to  the  sixth  power 
of  the  spatial  frequency.  Convolution  filters  with  large 
PSFs  (the  size  is  indicated  in  the  lower  left  of  each 
image)  and  the  frequency-domain  filter  give  better  point- 
target  enhancement  snd  rejection  of  the  edges  of  extended 
regions  of  clutter  (such  as  the  river). 
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FIGURE  16  -  The  high-pass  filtered  images  of  Fig.  15  after  application 
of  a  single-intensity  threshold  at  1^  +  5.5o  to  detect 
the  point  targets.  Convolution  filters  with  large  PSFs 
(the  size  is  indicated  in  the  lower  left  of  each  image), 
and  the  frequency-domain  filter  give  better  target-detec¬ 
tion  performance. 
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rates  of  2,  4,  S  and  16  per  image.  Figure  17  clearly  shows  that  a 
filter  with  a  large-area  PSF  (e.g.,  11x11  elements)  produces  a  better 
target-detection  capability  than  one  with  a  small-area  PSF  (e.g.,  3x3 
elements).  As  noted  in  Sect.  4.1,  in  the  present  case  the  target- 
detection  probability  is  a  nonlinear  function  of  the  number  of  detected 
targets,  so  the  observed  improvement  is  a  substantial  one.  For  the 
reason  described  later  in  this  section,  the  5x5-element  filter  gave 
poorer  results  than  the  3x3-element  one,  and  we  did  not  plot  them  in 
Fig.  17. 


As  the  images  of  Fig.  15  demonstrate,  one  reason  why  a  filter 
with  a  large  PSF  is  preferable  for  point-target  detection  is  because  it 
is  more  effective  at  attenuating  low-frequency  clutter  that  is  corre¬ 
lated  over  relatively  large  areas  of  the  image.  In  particular,  the 
edges  of  extended  objects  are  composed  mainly  of  high  spatial  fre¬ 
quencies  in  the  direction  perpendicular  to  the  edge,  but  considerable 
low  spatial-f requency  content  is  present  in  the  direction  of  the  edge 
Itself.  An  edge  can  be  removed  by  high-pass  filtering  only  if  the 
lowest  spatial  frequencies  that  compose  it  are  strongly  attenuated,  and 
this  requires  that  the  corresponding  dimension  of  the  filter  PSF  be 
comparable  to  the  effective  length  of  the  edge.  Intuitively,  it  is  the 
interleaved  positive  and  negative  filter  coefficients  evident  in 
Fig.  13  that  cause  the  cancellation  of  an  extended  edge. 

A  second  reason  why  a  filter  with  a  small  PSF  is  less  suitable 
for  point-target  detection  is  because  of  the  difficulty  one  has  in 
obtaining  a  frequency  response  that  is  circularly  symmetric.  Both  the 
point  targets  and  the  clutter  that  are  typically  present  in  satellite 
imagery  have  circularly  symmetric  spatial-frequency  distributions  cen¬ 
tered  on  the  origin.  The  high-pass  filter  required  for  clutter  rejec¬ 
tion  should  therefore  have  a  circularly  symmetric  frequency  response. 
However,  because  of  the  limited  number  of  coefficients,  and  the  rectan- 
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gular  grid  on  which  they  are  defined,  high-pass  filters  with  a  fre¬ 
quency  response  that  approximates  circular  symmetry  cannot  be  obtained 
by  convolving  a  small  PSF. 

This  Is  Illustrated  In  Fig.  18,  which  compares  the  frequency 
responses  of  the  five  present  convolution  filters  with  that  of  the 
desired  <o6  filter.  We  obtained  the  frequency  responses  by  applying 
each  filter  to  an  Image  containing  a  point  against  a  zero  background, 
and  then  calculating  the  modulus  of  Its  2-D  FFT.  The  gray  level  of 
each  display  Is  proportional  to  the  logarithm  of  the  frequency 
response,  with  full  white  representing  0  dB  and  full  black  -80  dB.  The 
zero  spatial  frequency  is  at  the  center  of  each  display. 

Note  the  prominent  lobe  structures  and  the  highly  directional 
nature  of  the  frequency  responses  of  the  filters  obtained  by  convolving 
small  PSFs.  Only  the  llxll-eleraent  filter  approximates  the  desired  u>6 
response.  The  presence  of  lobe  structures  and  other  deviations  from 
circular  symmetry  In  the  frequency  response  of  a  2-D  filter  will  Intro¬ 
duce  artificial  directional  structure  into  the  filtered  Images.  In 
particular,  the  asymetry  will  enhance  clutter  in  particular  spatial 
directions,  and  the  lobe  structure  will  enhance  clutter  that  has  spe¬ 
cific  periodicities. 

Examination  of  the  filtered  Images  of  Fig.  15  reveals  these 
artifacts.  The  5x5-,  7x7-  and  9x9-element  filters  have  produced  a 
preferential  enhancement  (l.e.,  a  reduced  attenuation)  of  clutter 
oriented  In  the  two  diagonal  directions  relative  to  that  with  a  verti¬ 
cal  or  a  horizontal  orientation.  In  addition  to  this  directionality, 
the  clutter  present  in  the  three  images  also  has  a  periodic  nature. 
This  Is  caused  by  the  local  response  maxima  and  minima  that  the  three 
filters  possess  In  the  mid-frequency  range.  In  comparison,  the  clutter 
present  in  the  Image  of  Fig.  15,  filtered  with  the  2-D  FFT,  shows  an 
Insignificant  amouit  of  directional  or  periodic  structure  and  the 
effects  are  only  slightly  evident  in  the  Image  filtered  with  the  11x11- 


FIGURE  18  -  The  frequency  responses  of  the  five  convolution  filters  and 
the  frequency -domain  one  employed  in  Fig.  15.  Zero  spatial 
frequency  is  at  the  center  of  each  image  and  the  two  axes 
correspond  to  the  horizontal  and  vertical  spatial  frequen¬ 
cies.  The  intensity  of  the  display  is  proportional  to  the 
logarithm  of  the  filter  response  with  full  black  corre¬ 
sponding  to  -60  dB.  The  small  area  convolution  filters 
deviate  significantly  from  the  desired  circular  symmetry. 
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element  PSF.  The  strong  presence  of  these  artifacts  In  the  Image  fil¬ 
tered  with  the  5x5-element  PSF  may  explain  why  the  point-target- 
detection  performance  of  this  filter  was  significantly  poorer  than  even 
that  of  the  3x3-element  filter. 


5.0  CONCLUSIONS 


Spatial  filtering  can  enhance  point  targets  located  in  a  clut¬ 
tered  background  in  images  obtained  from  a  satellite-based  mosaic  IR 
detector.  We  examined  by  simulation  and  compared  the  target-detection 
performances  of  a  variety  of  spatial  filters  implemented  by  direct 
convolution  and  by  2-D  FFT  filtering. 

When  a  single-intensity  threshold  is  used  to  detect  the  targets, 
it  is  preferable  to  use  a  spatial  filter  that  produces  considerably 
more  high-frequency  enhancement  than  a  conventional  Wiener  filter.  We 
found  that  a  filter  with  a  frequency  response  proportional  to  the  sixth 
power  of  the  spatial  frequency  was  effective  at  enhancing  the  point 
targets  and  rejecting  the  clutter  in  the  present  test  imagery.  For 
comparison,  the  corresponding  Wiener  filter  had  an  approximate 
frequency  response  between  the  square  and  the  cube  of  the  spatial 
frequency. 

Convolution  filters  with  large-area  PSFs  (e.g.,  11x11  elements), 
and  those  implemented  in  the  spatial-frequency  domain  with  a  2-D  FFT, 
are  more  effective  than  filters  with  small-area  PSFs  (e.g.,  3x3 
elements)  at  attenuating  the  edges  of  extended  clutter  regions.  In 
addition,  filters  with  small  PSFs  have  spatial-frequency  responses  that 
are  not  circularly  symmetric  and  which  can  possess  lobe  or  directional 
features.  These  properties  can  introduce  undesirable  artifacts  and  can 
reduce  the  clutter-rejection  effectiveness. 
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The  processing  required  to  detect  and  track  targets  from  a 
space-based  surveillance  system  will  include  temporal  and  multispectral 
processing,  as  well  as  the  spatial  filtering  of  Individual  image  frames 
considered  here.  In  future,  we  will  examine  how  temporal  and  multi¬ 
spectral  processing  can  be  incorporated  into  the  overall  processing. 
However,  such  processing  algorithms  will  be  particularly  dependent  on 
the  characteristics  of  the  imaging  system  (e.g.,  signal-to-noise  ratio, 
spatial  resolution,  frame  rate,  multispectral  response,  etc.)  and  on 
the  Imaging  scenario  (e.g.,  the  characteristics  of  the  target  and  the 
background,  viewing  geometry,  satellite  jit>0.r,  etc.).  Realistic 
imagery  obtained  with  a  space-based  mosaic  IR  detector  will  be  needed 
to  optimize  the  processing  required  to  perform  specific  target- 
detection  and  tracking  tasks. 
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laca,  etc.)  que  lea  flltrea  avec  une  petite  aurface  (e.g.  3  x  3  ou 
5x5  AlAaenta). 
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We  examined  by  simulation  the  oae  of  apatlal  filtering  for 
enhancing  the  vlalblllty  of  point  targeta  and  reducing  the  background 
clutter  In  Imagery  obtained  with  a  apace-baaed  staring  moaalc  I*  detec¬ 
tor  array.  An  experimental  digital  Image-proceaal ng  system  waa 
employed  to  Interactively  proceaa  Landsat  II  Images  containing  simu¬ 
lated  point  targeta.  First,  we  obtained  Improved  target  detection  and 
clutter  rejection  by  using  filters  that  produced  considerably  more 
high-frequency  anhanceawnt  than  the  conventional  Wiener  filter. 
Second,  convolution  filters  with  large-area  point-spread  functions 
(a.g.,  thoee  with  11  by  11  eleavnta  or  those  implemented  with  a  2- 
dlmens tonal  FFT  algorithm)  were  more  effective  st  attenuating  the  edges 
of  extended  regions  of  clutter  (e.g.,  those  corresponding  to  clouds, 
bodies  of  water,  etc.)  than  were  those  with  small  areas  (e.g..  3  by  3 
or  5  by  5  eleawnts). 
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We  examined  by  simulation  the  uae  of  spat  1st  filtering  for 
enhancing  the  visibility  of  point  targeta  and  reducing  the  background 
clutter  In  imagery  obtained  with  a  apace-based  atarlng  mosaic  IR  detec¬ 
tor  array.  An  experimental  digital  Image-proceaal ng  system  waa 
employed  to  interactively  process  Landast  II  liaagea  containing  simu¬ 
lated  point  targets.  First,  we  obtained  Improved  target  detection  and 
cluttrr  rejection  by  using  filters  that  produced  considerably  more 
high-frequency  enhanceawnt  than  the  conventional  Wiener  filter. 
Second,  convolution  filters  with  large-ares  point -spread  functions 
(e.g.,  those  with  11  by  11  eleawnts  or  thoee  Implemented  with  a  2- 
dimenslonal  FFT  algorithm)  weie  more  effective  at  attenuating  the  edges 
of  extended  regions  of  clutter  (e.g.,  those  corresponding  to  clouds, 
bodies  of  water,  etc.)  than  were  those  with  small  areaa  (e.g.,  3  by  3 
or  3  by  3  elements). 
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"Spatial  Filtering  for  Enhancing  Point  Targeta  In  Images  from  a  Spare- 
*aaed  Moaalc  IR  Detector" 
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We  examined  by  simulation  the  uae  of  apatlal  filtering  for 
enhancing  the  vlalblllty  of  point  targeta  and  reducing  the  background 
clutter  In  Imagery  obtained  with  a  apace-baaed  atarlng  moaalc  IR  detec¬ 
tor  array.  An  experimental  digital  imaRe-proceaalng  system  was 
employed  to  Interactively  proceaa  Landaat  11  Images  containing  simu¬ 
lated  point  targeta.  First,  we  obtained  Improved  target  detection  and 
clutter  rejection  by  using  filters  that  produced  considerably  more 
high-frequency  enhancement  than  the  conventional  Wiener  filter. 
Second,  convolution  filters  with  large-area  point-spread  functions 
(e.g.,  those  with  11  by  11  elements  or  those  implemented  with  a  2- 
dtmenslonal  FFT  algorithm)  were  more  effective  at  attenuating  the  edges 
of  extended  regions  of  clutter  (e.g.,  those  corresponding  to  clouds, 
bodies  of  water,  etc.)  than  were  thoae  with  small  areaa  (e.g.,  1  by  3 
or  3  by  3  elements). 
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Spatial  Filtering  for  Enhancing  Point  Targets  In  Images  from  a  Space- 
Baaed  Mosaic  IR  Detector" 
by  J.F.  Boulter 


We  examined  by  simulation  the  uae  of  spatial  filtering  for 
enhancing  the  visibility  of  point  targets  and  reducing  the  background 
clutter  In  Imagery  obtained  with  a  apace-baaed  atarlng  moaalc  IR  detec¬ 
tor  array.  An  experimental  digital  image-processing  system  waa 
e^>loyed  to  Interactively  process  Landaat  II  Images  containing  simu¬ 
lated  point  targets.  First,  we  obtained  Improved  target  detection  and 
clutter  rejection  by  using  filters  that  produced  considerably  more 
high-frequency  enhancement  than  the  conventional  Wiener  filter. 
Second,  convolution  filters  with  large-ares  point-spread  functions 
(e.g.,  those  with  11  by  11  elements  or  those  implemented  with  a  2- 
dlmenalonal  FFT  algorithm)  were  more  effective  at  attenuating  the  edges 
of  extended  regions  of  clutter  (e.g.,  thoae  corresponding  to  clouda, 
bodies  of  water,  etc.)  than  were  those  with  small  areas  (e.g.,  3  by  3 
or  5  by  5  elements). 


We  have  graphed  In  Fig.  17  the  number  of  targets  detected  In 
Landsat  Image  A  as  a  function  of  the  type  of  high-pass  filter  that  was 
used.  Results  are  shown  for  the  five  convolution  filters  considered  in 
Sect.  4.2  as  well  as  for  the  2-D  FFT  filter  employed  In  Sect.  4.1.  The 
graphs  were  obtained  as  described  in  Sect.  4.1  for  fixed  false-alarm 
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FIGURE  15  -  Landsat  Image  A  after  high-pass  filtering  with  five  dif¬ 
ferent  spatial-domain  convolution  filters  and  a  frequency- 
domain  one  with  a  response  proportional  to  the  sixth  power 
of  the  spatial  frequency.  Convolution  filters  with  large 
PSFs  (the  size  is  Indicated  In  the  lower  left  of  each 
Image)  and  the  frequency-domain  filter  give  better  point- 
target  enhancement  and  rejection  of  the  edges  of  extended 
regions  of  clutter  (such  as  the  river). 
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FIGURE  16  -  The  hlgh-psss  filtered  images  of  Fig.  15  sfter  application 
of  s  single-intensity  threshold  st  1^  +  5.5o  to  detect 
the  point  targets.  Convolution  filters  with  lsrge  PSFs 
(the  size  is  indicated  in  the  lower  left  of  each  image), 
and  the  frequency-domain  filter  give  better  target-detec¬ 
tion  performance. 
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